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ABSTRACT 

We analyse Herschel/SPIRE images of the edge-on spiral galaxy NGC 891 at 250, 350, and 500 /urn. Using a 3D radiative transfer 
model we confirm that the dust has radial fall-off similar to the stellar disk. The dust disk shows a break at about 12 kpc from the 
centre, where the profile becomes steeper. Beyond this break, emission can be traced up to 90% of the optical disk on the NE side. 
To the SW, we confirm dust emission associated with the extended, asymmetric HI disk, previously detected by the Infrared Space 
Observatory (ISO). This emission is marginally consistent with the large diffuse dust disk inferred from radiative transfer fits to optical 
images. No excess emission is found above the plane beyond that of the thin, unresolved, disk. 
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1. Introduction 

The dust distribution in the edge-on galaxy NGC 891 has been 
studied in detail. Radiative transfer fits to optical/near-infrared 
(NIR) images have shown that its extinction lane can be de- 
scribed with a dust disk that is radially wider but vertically thin- 
ner than the stellar disk. (Because its exponential radial scale- 
length is 1.5-2 times and the vertical scalelength half that of the 
stars, so her eafter we call it the diffuse dust disk; Xilou ris et al.l 
I1998L 119991 -1 This disk is found to have moderate extinction 
properties, since the central optical depth perpendicular to the 
disk is about 1 in the B-band. However, this is not large enough 
to explain the amount of energy emitted by dust and observed 
in the far-infrared (FIR) and submm spectral energy distribu- 
tion (SED): the energy balance requires that up to three times 
more starlight is absorbed by dust than predicted by the dif- 
fuse dust disk ( Po pescu et al.l 120001) . Similar dust disk proper- 
ties and the same energy balance problem have been found 
in other edge-on galaxies jXi louris etaD[l999; MisiriotisetaD 
l200UlDasvra et alj|2005t lBianchill2007L l2008fc iBaes et alj|2010h . 

Radiative transfer models of dust emission in NGC 891 solve 
the energy balance problem by including a second dust distribu- 
tion associated with the molecular gas, of mass comparable to 
that of the diffuse dust disk. This second component escapes de- 
tection in optical images b ecause it is thinner than diffuse dust 
(Pope scu et al.ll200fl 12.01 ll) or it is clumpy in nature, thus pro- 
ducing effects that are less coherent than the extinction lane 
caused by diffuse dust (Bianchi 2008). Hereafter, we call this 
second compone nt the clumpy d ust disk. A model for NGC 891 
was presented in lBianchil (120081) . made with the radiative trans- 
fer code TRADING. Though the model includes several fea- 
tures, we found that emission beyond 100/vm is caused mainly 
by an exponential disk of evolved stars (radial scalelength 4 
kpc and vertical sca lelength 0.4 kpc, as derived in the NIR by 
IXilouris et al.l ll999) heating at thermal equilibrium both the dif- 



fuse dust disk (radial scalelength 8 kpc and vertical scalelength 
0.2 kpc) and the external envelope of the clouds in the clumpy 
dust disk (with clouds distributed exponentially with the radial 
scalelength of the stars, and the vertical scalelength of the dif- 
fuse dust disk). 

In this Letter we compare the structural properties of the 
dust distribution in NGC 891 inferred indirectly from extinc- 
tion studies and/or model predictions with those deriv ed from re 



cent s ubmm images taken by the SPIRE instrument ( Griffin et al 



2010) aboard the Herschel Space Observatory (Pilbrattetal 



2010). In particular, the high sensitivity of SPIRE, together with 



the edge-on configuration of the galaxy, allow studying submm 
emission up to large radial distances from the galactic centre. 
This is needed to directly verify the presence of the diffuse dust 
disk, whose predicted emission could become visible over that 
of the clumpy dust disk only in the outskirts of the galaxy. 



2. Observations and data reduction 

SPIRE photometric observations at 250, 350, and 500 /itn were 
obtained as part of the Guaranteed Time Key Project Very 
Nearby Galaxy Survey (P.I. C. Wilson). The galaxy was ob- 
served in large map mode, covering an area of 20'x20' centred 
on the object with two crossscans, using a 30" s _1 scan rate. 
Data was retrieved from the Herschel S cience Arc hive and re- 
duced with the dedicate software HIPB3 (Ottl l2010l) . As the cur- 
rent pipeline simply removes median baselines from the time- 
lines, causing artifacts in the maps from the presence of bright 
extended sources, we subtracted from each bolometer timeline a 
linear fit obtained after carefully masking the data over the re- 
gion covered by the galaxy. Maps were then produced using the 
naive mapmaking procedure within HIPE. The resulting images 
(Fig- LU nave uniform backgrounds. We also experimented with 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



1 We used HIPE 5.2 and the SPIRE calibration tree v. 5.1. SPIRE 
characteristics (point spread function - PSF - FWHM width and beam 
area, empirical PSF's, calibration uncertainty, and colour corrections) 
are taken from the SPIRE Observers' Manual (v2.3, 201 1). 
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Fig.l. SPIRE images of NGC 891 at 250 

(left), 350 (centre), 500yum (right). Each panel 
shows an area of 14'x28' centred on the 
galaxy's radio continuum coordinates a = 
2 h 22 m 33 s .O, 6 = 42° 20' 57/2 (J2000.0; 
lOosterloo et al. 2007). The main beam size is 
indicated by the filled circle (FWHM: 18'.'1, 
24'.'9, and 36'.' 4). Images have pixel sizes of 6, 
8, 12" at 250, 350, and 500 yum, respectively 
(i.e. about 1/3 of the measured FWHM). The 
crossed lines in the 500/jm image show the 
scan directions. The sky rms noise is a = 0.74, 
0.44, and 0.23 MJy sr -1 , from 250 to 500 yum. 
Contours are shown at 3, 10, 100, 200-cr over 
all images, at 500-cr in the 250 and 350yum 
map, at 1000-<r only in the 250/um map. North 
is up, east to the left. 



median baseline remov al and the alte rnative mapmaking soft- 
ware Scanamorphos v7 dRousse but found no significant 
difference for any of the analyses presented here. 

3. The data vs the model 

Total flux densities are 169+25, 72+11, and 26±4 Jy, at 250, 
350, and 500 fan, respectively. These values include a colour 
correction appropriate to the observed spectral index (F v oc v 2 ' 7 ): 
since the source is only resolved along the major axis (see later), 
we used a correction halfway between that for point and ex- 
tended sources. (We multiplied the pipeline flux densities by 
0.94, 0.95, 0.94.) We adopted a very conservative 15% cali- 
bration uncertainty for extended emission (see the Observers' 
Manual, Sect. 5.2.13). Other sources of uncertainty (photome- 
try, point source/extended calibration) are much smaller, ^ 4%. 

SPIRE flux densities are shown in Fig. [2] together with the 
FIR data at 60 and lOOyurn from the IRAS satellite and a t 170 and 
200/um from the ISO satellite (see iPopescu et alJl20lH and ref- 
erences therein), as well as the submm flux densities at 350, 550, 
and 850/mi from the HFI instrument aboard the Planck satellite^. 

Major axis profiles are shown in Fig. [3] To increase the 
signal-to-noise ratio, they were averaged over two beams across 
the plane. Residual sky gradients were estimated from two par- 
allel strips on both sides of the galactic plane, well beyond the 



2 Since NGC 891 is marginally resolved by Planck, we used 
the flux densities determined by fitting the data with a 2-D 
Gaussian m odel, as reported in the Early R elease Compact Source 
Catalogue dPlanck C ollaboratio n et al.l 1201 lh . The quoted calibra- 
tion uncertainties are 7% at 350 and 550/zm and 2% at 850 jjm 
dPlanck HFI Core Team et aDl201 lh . Colour corrections for a y 2 7 spec- 
trum were taken from the Catalogue Explanatory Supplement. We did 
not correct the flux density at 850yum for the contamination due to the 
CO(3-2) line. Th is was found to be about 5% in SCUBA observations 
jlsraeletal.ll 19991) . and it is likely to be lower for the wider 353GHz 
bandpass. The Planck flux density at 350/jm appears to be higher than 
SPIRE, though still within 1-cr, for the large SPIRE calibration uncer- 
tainty adopted here. We do no t use other submm observations fro m the 
literature (for references, see: Pocescu et al. 2011; Bianchi 2008), be- 
cause they could have been maimed by a limited spatial coverage and 
uncertain PSF sidelobe corrections. We note, however, that the Planck 
flux density at 850/um is within th e errors of SCUBA observations 
dAltonetal.ll 19981 : llsrael etal.lll999h . 
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Fig. 2. The integrated FIR-submm SED of NGC 891. 



emission perpendicular to the major axis. They were found to 
be smaller than 0.5-cr in all bands, so they do not affect the pro- 
files significantly. Profiles were extracted using a position an- 
gle of 22.9°, w hich was determined with the fitting procedure of 
Bianchi (2007). The inner profile of the galaxy can be roughly 
described as the edge-on projection of an exponential disk, with 
radial scalelength increasing from 4.6 kpc at 250,um, to 5.1 kpc 
at 350/mi, a nd to 5.2 kpc at 500/um. (W e assumed a distance 
D=9.5 Mpc: lvan derKruit&Searldll981l )A break occurs in all 
bands at about 12 kpc from the centre (4 '4), after which the pro- 
files become steeper (with a radial scalelength of about 1 kpc). 
Such radial breaks are analogo us to what is found in most stellar 
disks dPohlen & Tru iillo 2006) and have already been observed 
with SPIRE dPohlen et alj|2010l) . Though the change in slope in 
NGC 891 is quoted as appearing at about 6' (co ming to a sharp 
truncation at 7'5: Ivan der Kruit & Searlell 19811) . the optical im- 



ages used in Xilouris et al. (1999) indicate a break similar to 



the one detected for dust. SPIRE clearly improves over previous 
submm observations. SCUBA was able to d etect the dust emis - 
sion at 850//m within the radial break only (lAlton et al.l fl998). 
The 850yum morphology for the inner region is very similar to 
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Fig. 3. Surface brightness profiles along the 
major axis from SPIRE images (solid lines) and 
the corresponding models by Bianchij d2008h 
(dot-dashed lines) at 250 (top, blue, lines), 350 
(middle, green, lines), and 500pm (bottom, red, 
lines). We also show the profile for the model 
that includes dust heating from the EBL (dot- 
ted line, for 250pm only). For each band, pro- 
files are averaged over two beams perpendicu- 
lar to the major axis. They are shown for flux 
densities > 1.2, 0.75, and 0.45 MJy sr~' for 
250, 350, and 500pm, respectively, correspond- 
ing to S/N > 3 for the average. Negative dis- 
tances from the centre refer to the SW end 
of the galax y. The optical size of the galaxy 
(D 25 = 13!5: lde Vaucouleurs et aTlll991l) is in- 
dicated. See text for other details. 



SPIRE ima ges, and a radial scalelength of 5.3 kpc was fitted to 
the profile (lAlton et al.l l2000). close to what we find at 500pm. 

After the break, the profile on the NE side of the galaxy 
can only be traced up to 6', or 90% of the optical size (as- 
suming D 25 = 13f5 = 37 kpc; |de Vaucouleurs et al .1 1199 lh . On 
the SW end, instead, significant emission is seen in excess of 
the break, up to the optical size and beyond. The excess has 
already been detected i n ISOPHOT images at 170 and 200pm 
(Porjescu & Tuffs 2003); it corresponds to the SW extension of 
the asymmetric HI disk, as can be seen by comparing the 3-<x 
SW appendage of Fig. [T] with the 10 21 cm" 2 HI column density 
contou r at a comparable resolution in Fig. 1 of lOosterloo et alJ 
(l2007h . From the data, we derived a dust mass for this feature of 
6 x 10 s M© (though with an uncertainty of almost a factor two, 
mainly because of the difficulty deriving temperatures using only 
SPIRE fluxes). The HI mass corresponding to th e SPIRE detec- 
tion ( estimated from the column density levels in Oos terloo et alJ 
l2007t) is about 10 8 Mo, resulting in a dust-to-gas mass ratio of 
0.006, similar to what is found for the i nner disks of other galax- 
ies and of the Milky Way (see, e.g., iDraine etal 1 120071) . This 
confirms that the gas is not pristine, and it supports the view that 
the lopsidedness of NGC 89 1 originated from the perturbation of 
the ga lactic disk during a fly-by of the nearby companion UGC 
1807 dMapelli et al.ll2008l) . The HI mass for the whole SW disk 
extension is 2.5 x 10 s M©. Assuming the same dust-to-gas mass 
ratio, we estimated a total dust mass of 1.5 x 10 6 M© for t he 
region. (A similar result was found bv lPopescu& Tuffs 200l) 

Previous radiative transfer models of NGC 891 underpredict 
the new s ubmm SEP. Fo r instance, the 250pm flux density pre- 
dicted bv lBianchl (120081) is about 45% lower than observed (see 
Fig. [2] where we also s how the similar, but independent, model 
of iPopescu et all 1201 lh . We revised the model to describe the 
submm SED and the radial profiles observe d by SPIRE better. 
Following the procedure of Bianchi (2008), we increased the 
mass in the clumpy dust disk (from 0.5 to l.lxlO 8 M©, making 
up 70% of the total dust mass) and fine-tuned the other parame- 
ters (mainly, we raised the galaxy intrinsic luminosity from 7.8 
to 8.5 x 10 s LQ) to reproduce the stellar SED as well. The SED of 
the new model is shown in Fig. [2] It is only marginally consistent 
with the 100pm IRAS flux density, and it does not reproduce the 
170 and 200pm flux densities from ISO. With the FIR/submm 



emissivity of the adopted dust grain model (IDraine "&Lill2007l) . 
and for the interstellar radiation fields implied by the relative ge- 
ometries of dust and stars assumed within the radiative transfer 
code TRADING, it was not possible to find a model able to de- 
scribe all datapoints. 

Modelled major axis profiles (Fig. [3} were obtained by con- 
volving simulated images in the three bands with the respective 
empirical PSFs, regridding them on the chosen pixel size and 
averaging the profile as for real data. For a match of the broad 
surface brightness distribution along the major axis (as well as 
of the SED), we had to use the same radial scalelength (5.7 kpc) 
both for the clumpy dust disk and for the stellar disk. The sec- 
ond result is puzzling, since most of the observed stellar emis- 
sion (and a good fraction of the absorbed energy in the current 
model) comes from the NIR, for which the radiati ve transfer fits 
found a smaller stellar radial scal elength (4 kp c; IXilouris et alJ 
fl999l) . Also, as already noted in Bianchi (2008), emission from 
the clumpy dust disk appears smoother than the real data, reveal- 
ing that dust inhomogeneities on a larger scale than molecular 
clouds (i.e. rings, spiral arms) should be considered in models. 

The break in the profile is reproduced in the model by a 
sharp cut of both the clumpy dust disk and the stellar disk at 
13 kpc from the centre (4.7' in Fig.[3]B Beyond the break, the 
modelled emission is only due to the diffuse dust disk derived 
from optical observations. The mean dust temperature in this 
disk is low (10 < T < 13K) because the heating sources (mainly 
stars at the edge of the disk; the bulge contribution is negligible) 
are no longer mixed with dust. The predicted surface brightness 
matches the observations only on the SW, where emission comes 
from the excess prev iously described. The axisymmetric diffuse 
dust disk derived by Xilouris et al.l (1 19991) could thus mimic the 
contribution to extinction of the (asymmetric) extended HI com- 
ponent. However, it is to be noted that the model predicts a lower 
temperature than what is derived in the SW appendage (15K), 
and the mass of the diffuse dust disk beyond 1 3 kpc (8 x 1 6 M©) 
is 5x higher than what is estimated from the HI data, which is 
too high a factor to be explained by the uncertainties in the tem- 
perature and in the various model assumptions. Also, model pre- 



3 Conversely, the break could be produced by truncating the clumpy 
dust disk only, if the model did not include a diffuse dust disk. 
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Fig. 4. Data and (scaled) model profiles at 250ji/m perpendicular 
to the major axis through the centre of the galaxy and the two 
secondary peaks. The profile through the SW peak has been di- 
vided by 10, that in the NE by 100). Negative heights refer to the 
NW side of the disk. 



dictions are likely lower limits to the real emission from a diffuse 
dust disk. If stars extend more than what is needed to fit the radial 
break, the emission will be higher. Instead, if grains are isolated 
from stars, the extragalactic background light (EBL) might con- 
tribute significantly. The dotted line in Fig. [3] shows the effect 
on the 250jt/m profile of the EBL PI .0 model of lAharonian et al.1 
(2006). Emission from the diffuse dust disk in this region of the 
model could rise by a factor twcQ. 

Finally, Fig.|4]shows the 250^/m profiles perpendicular to the 
major axis, for three cuts passing through the galaxy centre and 
the two secondary peaks at about V from the centre on the NE 
and SW. Data were averaged over two beams perpendicular to 
the cut, and a residual background was subtracted from regions 
beyond the emission. Model profiles were scaled to match the 
data on the major axis. Emission above the plane is almost en- 
tirely dominated by the PSF Airy rings picking up light from 
the major axis (because the adopted dust vertical scalelength, 
0.2 kpc or 4", is much smaller than the main beam). Results are 
similar in the other SPIRE bands. Excess emission is only seen 
on the NW side of the minor axis profile, in coincidence with 
the background X-ray source CXO J022224.4+422138, which 
contributes to most of the excess; however, part of it could come 
from another background source aligned with the minor axis or 
from a filamentary structure protruding from the galactic cen- 
tre. Excluding this feature, SPIRE images do not show evidence 
of any diffuse halo emi ssion. This is at odds with the previous 
Spitzer observations of Bur gdorf et al.l (120071) . which detected 
MIR emission extending up to 5 kpc from the plane. At 22//m, 
a surface brightness of 0.1-0.2 MJy/sr is found at 1-1 'A along 
the minor axis. Instead, at 250jum, residual emission from the 
halo is less than 1 MJy/sr. Such MIR/submm flux density ratios 
can be produced by typical Milky Way dust, if the interstellar 



4 The temperature of dust heate d by the EBL only is about 1 1 K, 
adopting mean dust properties from Draine & Li (2007|). 



radiation field is l Ox higher than the local (see, e.g., Fig. 13 in 
Drai ne &"l1|2007I) . These heating conditions might be compat- 
ible with the hypothesis that the detec ted MIR radiation come 
from dust surrounding halo AGB stars (Burgdorf et al. 2007). 

4. Summary & conclusions 

We have analysed SPIRE observations of the edge-on spiral 
galaxy NGC 891 and examined both their contribution to the 
galaxy's SED and the galaxy's morphological properties at these 
wavelengths. Emission within the central 4'4 broadly follows 
a radial exponential distribution with a scalelength similar to 
those derived for th e stars in the surface brightness fits of 
IXilouris et al.1 (1 19991) . Emission above the galactic plane can be 
explained with radiation from the galactic plane picked up by the 
PSF wings. The distribution of dust along the vertical direction 
is thus thinner than the SPIR E PSFs, and compa tible with the 
vertical scalelength derived bv lXilouris et al.l dl 999b . 

Beyond A'A from the centre, the radial profile shows a break 
and becomes steeper. On the NE side, the submm disk comes 
to an end within the optical size of the galaxy. Instead, on the 
SW side the profile decline is interrupted at about the optical ra- 
dius, where significant emission is seen in all bands. With their 
high resolution, SPIRE observations confirm the excess emis- 
sion associated with the asymme tric HI disk, already f ound in 
deep ISOPHOT observations bv lPopescu & Tuffs (2003|). 

The emission within the central 5' can be attributed to the 
clumpy dust disk required by radiati ve transfer mo dels to fit 
the FIR/submm SED of the galaxy dBianchil 12008). The HI- 
associated dust emission in the SW instead has surface bright- 
ness levels that are marginally consistent with those predicted 
for the diffuse dust disk, an important dust component (of simi- 
lar mass to that of the clump y dust disk) respo nsible for the op- 
tical/NIR extinction lane dXilouris et al.lll999f) . Since extended 
diffuse dust disks are routinely found in optical/NIR radiative 
transfer fits of edge-on galaxies, we hope that further Herschel 
observations of these objects (as those planned within the ap- 
proved projects HEROES and NHEMESES), together with ra- 
diative transfer models and deep observations of their atomic gas 
disks, will help for understanding if the diffuse dust disk is in- 
deed a component associated with the atomic gas and is a feature 
common to all spiral galaxies. 
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